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Abstract: In this paper we study the nonperturbative structure of the axial current 
and evaluate the two-point distribution amplitudes f e~ tq '^ (0\ti;(x)T^(y) J^(£)|0) in the 
framework of the instanton vacuum model in the leading order in O (N c ). We perform a 
direct numerical test of the relations between the axial current and the pion distribution 
amplitudes, imposed by PCAC, and found excellent agreement. 
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1 Introduction 

According to the Standard Model, the interaction of W- and Z-bosons with quarks has a 
V — A structure. While the hadronic structure of the vector current has been well estab- 
lished in photon-hadron and photon-nuclear processes, the structure of the axial current 
is less known. Due to the spontaneous breaking of the chiral symmetry the structure of 
the axial current is different from the vector case. In order to describe the longitudinal 
component of the axial current a widely used tool is a phenomenological PCAC relation, 
which was proposed in the pre-QCD era [1-4], and relates the longitudinal part of the 
interaction amplitude of the axial current with that of the pion [5]. However, PCAC is not 
applicable to a description of the transversely polarized axial current. 

The standard description of high-energy processes is based on factorization [6] of the 
corresponding amplitude to a process-dependent hard part, evaluated in pQCD, and a 
target-dependent soft distribution amplitude (DA), which is either extracted from fits to 
data, or evaluated relying on low-energy models for the target. 

In this paper we evaluate the two-point quark distribution amplitudes (DAs) of the 
axial current in the framework of the instanton vacuum model (see [7-9] and references 
therein). The important advantage of this model is that it has a built-in dynamically broken 
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chiral symmetry and as a consequence, reproduces the low-energy chiral structure of QCD. 
Recently this model was used for evaluation of the vector current hadronic structure [10, 11], 
pion distribution amplitudes [12-17], correlators of vector and axial currents [10, 18, 19], as 
well as different low-energy constants. The nonperturbative structure of the axial current is 
particularly important for processes which include soft kinematics or large distance between 
quarks. 

The paper is organized as follows. In Section 2 we briefly overview the basic elements 
of the instanton vacuum model (IVM), which is used for further evaluations. In Section 3 
we give definitions for the DAs, and in Sections 3.1-3.8 provide analytic expressions and 
numerical results for each of the DAs. In Section 4 we directly check the accuracy of the 
PC AC relation within the IVM model. We summarize observations and make conclusions 
in Section 5. 



2 Instanton vacuum model 

The central object of the model is the effective action for the light quarks in the instanton 
vacuum, which in the leading order in iV c has the form [8, 9] 

S = J d 4 x (^y In A + 2<$> 2 {x)i) (p + v + a-ys-m - cLf ® $ • T m ® fL) tjj) (2.1) 

where T m is one of the matrices, T m = 1, if, 75, or if 75; tp and <!> are the fields of constituent 
quarks and mesons respectively; N/V is the density of the instanton gas; v = v «7^ is the 
external vector current corresponding to the photon; L is the gauge factor, 

L(x,z) =Pexp(i j ddv^O + a^C)^ , (2.2) 

L(x,z) = 7 L(x,z) t 7o (2.3) 

which provides the gauge covariance of the action [20, 21]; and f(p) is the Fourier transform 
of the zero-mode profile in the single-instanton background. In this paper we employ the 
dipole-form [8] for the parameterization of the formfactor, 

f<P) = (2-4) 

with L ~ 850 MeV. It is worth mentioning that the distribution amplitudes, due to their 
definition as nonlocal quark operators with light-cone separation, are much more sensitive 
to the choice of the formfactor compared with various vacuum condensates. Varying the 
formfactor f(p) completely different results were obtained for the leading- twist pion DA 
in [12, 13, 22]. 

In the leading order in N c , we have the same Feynman rules as in perturbative theory, 
but with a momentum-dependent quark mass fj,(p) in the quark propagator 

S(p) = — . (2.5) 
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The running mass of the constituent quark has a form 



Hip) =m + Mf(p), (2.6) 

where m ~ 5 MeV is the current quark mass, M ~ 350 MeV is the dynamical mass 
generated by the interaction with the instanton vacuum background. Due to presence of 
instantons the vector current - quark coupling is also modified, 

ij = ty/* ^V = v + V nonl , (2.7) 
a = a n ^^A = a + A nonl , (2.8) 

In addition to the vertices present in perturbative QCD, the model has nonlocal terms with 
higher-order couplings of currents and mesons. The exact expressions for the nonlocal terms 
V nonl ,A nonl depend on the choice of the path in (2.2), and one can find in the literature 
different results [10, 11, 13, 23]. This ambiguity in the DAs is important for the transversely 
polarized, but is negligibly small for the longitudinally polarized axial current. In what 
follows we use the parameterizations, 

Vnonl = Vn (* M ^f (/ faf ~ f fo) 2 )) , (2-9) 

A nonl = (iM^±^ (/ (pi) - / (p 2 )) 2 ) , (2.10) 
V Pi P2 / 

where pi, p 2 are the momenta of the initial and final quarks. 

3 Distribution amplitudes for the axial current 

The DAs of the axial current are defined via 3-point correlators, 

/ Ae-^(0|^(y)r^(x)j5(O|0>, (3.1) 



where x and y are light-cone coordinates of the quark and antiquark, q is the momentum 
flowing through the axial current and T is one of the Dirac matrices, as was defined in 
(2.1). The structure of the axial current is different from the vector one because of the 
spontaneous chiral symmetry breaking and existence of near-massless pions. In particular, 
the axial current can fluctuate into a pion prior to the production of a qq pair. Therefore, 
the correlator (3.1) consists of two terms, schematically presented in the Figure 1. 

One term comes from the combined contribution of the intermediate heavy states (oi 
meson, 3ir, etc.), and the other one represents fluctuations of the axial current into a pion. 
The chiral symmetry embedded into the model relates the two terms as, 

^ = y(bulk) + $ (p»n) = f _ ^ Ulk \ (3.2) 

This form of the DA explicitly satisfies PC AC. In what follows we concentrate on the part 
of the amplitude presented in the dispersion relation for the amplitude by the bulk of heavy 
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Figure 1. The DA contains two terms corresponding to either intermediate heavy axial states 
(left), or to a pion (right), which we label by bulk or pion respectively. 

states excluding the pion pole (left pane of Figure 1) [24, 25], tacitly assuming that the full 
DAs are given by Eq. (3.2). We define these DAs as, 
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where q is the 4- momentum carried by the axial current; u and u = 1 — u are the fractional 
light-cone momenta carried by the quark and antiquark; z = x — y; p^ is the "positive 
direction" vector on the light-cone, is the "negative direction" vector on the light- 



cone, with normalization conditions p ■ n = 1, p 



n 



0. Transverse dimensions are 



chosen in such a way that the vector q has only components in the p, n- plane. Without 
any loss of generality, in what follows we choose a system where g+ = q ■ n = 1. The 
normalization constant Ja is a dimensional parameter introduced in order to make the 
distribution amplitudes dimensionless. Its value is fixed by the condition 



f da $|| (a,Q 2 = 0) = 1. 
Jo 



(3.7) 
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If we define an "effective" axial meson state vector |^4( A )(g)) as 
A^(q)) = [ ^xe^- x ef\q)Jp(x)\Q), 



where the polarization vectors are defined as [26] 



(3.8) 



(3.9) 



n ■ e C^0 = 0, 



(<?) 



■1, 



(3.10) 
(3.11) 



then we may rewrite Eqns. (3.12-3.15) in a standard form as distribution amplitudes of the 
effective axial meson state [27], 



(0\$(V) W>P(x)\A(q)) = if A / dae^- a > z x 

J o 



e< A > • z 
(p ■ z) 2 



93(a) 3,-12) 



<0|^(y) 7 ^(^)|A(g)) 



^ Pp zj\ae^-^ 9 -^P- (3. 
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(0 \i> (y) (x)\ A(q)) = f A jf da e^-")** (( e (A=^ _ e ( A =^) $ ± (a) + 

(3.14) 
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p ■ z 



■h 3 (a) , 



(0${y) l5 iP(x)\A(q)) = f A e W -n f daS°^ a > z ■ 

Jo 



h<*\ a ) 



(3.15) 



The distribution amplitudes $||(n), &±(u) are of twist-2; <?j° , <7j" , /iff , /iff are of twist- 
3; 53, ^3 are of twist-4. Chiral parity: all wave functions in (3.12), (3.13) are chiral even, 
all wave functions in (3.14), (3.15) are chiral odd. 

Next step is a numerical evaluation of different DAs. 
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3.1 Evaluation of $n (a, q 2 



In order to extract this function, we have to take a convolution of (3.12) with ppn^, which 
yields 



8N C f drd 2 l ± 



-|) 7+75^ (|) Ji(x) 



fA J (2tt) 



H{l)n{l + q) + l\ + aag 2 

.a 2 +M 2 (0) (a+9) 2 +/i 2 (/+(?)) 

M (/(/ + g) - /(Q) 2 (2/~ + g 2 q) (/i(Z)q + jiil + g)q) 



(3.17) 



Obviously, this function is symmetric under the a — > 1 — a transformation. The 
normalzation condition (3.7) gives for /a 



fA = 8N C 



d 4 l ± 
(2kY 



M{f{l)fl 



(i 2 + fi 2 (i) f i 2 (i 2 + fi 2 (i)y 



As is discussed in the Section 4, if the PC AC relation is valid, then Ja ~ V^f 2 . 
In Figure 2 the DA $n(a) is plotted as function of a for several values of Q 2 . 
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Figure 2. Distribution amplitude $11(0) vs. a for several values of Q 2 



3.2 Evaluation of $j_ (a, q 2 ) 

In this section we evaluate the function 
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/ 2 
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(3.19) 



(Z + <?) 2 - Z 2 ( Z 2 + ^(J)) ( (/ + ?) 2 + ^ + ?) 



l+=- aq -\ 



This function is antisymmetric under the a — > 1 — a transformation. In Figure 3 the 
distribution amplitude is shown for several values of Q 2 . 




Figure 3. Distribution amplitude $j_(a) for several values of Q 2 . 



3.3 Evaluation of gy (a, q 2 ) 

In ordi 
yields 



In order to extract this function, we have to take a convolution of (3.12) with \g^^ which 



g± q ) 



8N C fdl-d 2 l± 



if A J (2irY 



+ q) + (2a - l)l-q+ + a* 
(Z 2 + /, 2 (Z)) ({l + qf+^d + q) 



(3.20) 
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72 



M{f(l + q)-f{l)f ( M (Z + g)- M (0) 



(I + q) 2 - I 2 ( |2 + M 2( Z )) ^ + q f + ^ + q) 



This function is symmetric under the a — > 1 — a transformation. In Figure 4 the 
distribution amplitude is shown for several values of Q 2 . 
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Figure 4. Distribution amplitude <jr" (a) for several values of Q 2 . 
3.4 Evaluation of gy (a, q 2 ) 

In order to extract g^ (u), we have to take a convolution of (3.13) with e^pvPp'Wu' , which 
yields 



<#> (a, g 2 ) = - 2£ "y g | dSe"** | ^0 



16iY r fdl~d 2 h 



fA J (2vr)4 



(3.21) 



\ J A l+=—aq+ 

This function is antisymmetric under the a — > 1 — a transformation. In Figure 5 the 
distribution amplitude is shown for several values of Q 2 . 

3.5 Evaluation of hj\' (a, q 2 ) 

In this section we evaluate the function 



h® (a, q 2 ) 



{VvUp - p p n y ) P/3 
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((I + g )a-Za)(P + ^(Z))((I + (? ) 2 + Ai a (H-«); 

\ / J (+=— aq+ 

This function is antisymmetric under the a — > 1 — a transformation. In Figure 6 the 
distribution amplitude is shown for several values of Q 2 . 
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Figure 6. Distribution amplitude h,, (a) for several values of Q 2 . 
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3.6 Evaluation of h!^ (a, q 2 ) 

In order to extract (u), we have to take a convolution of (3.15) with 2pp/ /a, which 
yields 



hf (a, q 2 ) 



' fA 



I ^ e i(a-0.5)p-z t 







WN r [ dl~d 2 li 
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+ 



l]_ + {\-2a)l -a^ + KW + Q)) M (f (I + q) - f (I)) 2 (2l~ + q 2 a) 



{l + q) 2 -l 2 



This function is symmetric under the a — > 1 — a transformation. In Figure 7 the 
distribution amplitude is shown for several values of Q 2 . 

i — i 1 1 1 1 1 1 1 1 — t 

Q 2 = GeV 2 
Q 2 = 0.5 GeV 2 
Q 2 = 1 GeV 2 
Q 2 = 2 GeV 2 
Q 2 = 4 GeV 2 



2.5 



. 2.0 



©) 

B 



7.5 



-5 



1,0 



0.5 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
a 



Figure 7. Distribution amplitude nj? (a) for several values of Q 2 . 



3.7 Evaluation of g% (a, q 2 ) 

In order to extract this function, we have to take a convolution of (3.12) with ppPfj,, which 
yields 
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This function is symmetric under the a — > 1 — a transformation. In Figure 8 the distribution 
amplitude is shown for several values of Q 2 . 
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Figure 8. Distribution amplitude 33(a) for several values of Q 2 



3.8 Evaluation of /13 (a, q 2 ) 

In this section we evaluate the function 



(«, q 2 ) = ^( 9 ^P»-gp»P») j d'xe-^j g e ^-0-%-( |^(-|n 
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2 {I + qf - P (/2 + ^ h x + q f + M 2( Z + ? ) N 

V / J / + = — ciq^ 

This function is antisymmetric under the a — > 1 — a transformation. In Figure 9 the 
distribution amplitude is shown for several values of Q 2 . 
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Figure 9. Distribution amplitude /13(a) for several values of Q 2 . 
4 Tests of PCAC 

The PCAC hypothesis was proposed in the pre-QCD era [1-4] and has been intensively 
used as a phenomenological tool for describing the longitudinal part of the axial hadronic 
current at small virtualities. In operator form, the PCAC relation is 

d^ a = Umlr, (4-1) 

where J^' a is the axial current; (j) a is the effective pion field; and /„- is the pion decay 
constant. 

The pion DAs are defined as [25, 28], 

(0|$(y)7M->^(s)H?)> = if*V2 [ dae^-y+w^ x 

J 



{V\i>{y)l^(x)\v{q)) = -if v V2 m } I* dae^y+^U'glia), (4.3) 

m u + m d J 



2 ri 

(0 \i> (y) (x)\ n(q)) = ~f„y/2 ^ / dae i ^+^ x 

3 m u + m d J 

x (jppZv-pvzJ^lia). (4.4) 

p ■ z 

These expressions can be bridged via Eq. (4.1) with the corresponding DAs of the axial 
current, Eqs. (3.17), (3.22), (3.23) and (3.24), resulting in, 
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f A $\\ (a, q 2 = ml) = fiV^fo.^a) (4.5) 

(a, <7 2 = "4) = - JL ^—ip4;ir (a) (4.6) 

/a/i | ; (a, 9 = = -— - <^(a) (4.7) 

f A hf (a, q 2 = ml) = ^?g£$> ( a ) (4.8) 



The DAs of the pion Eqs. (4.2)-(4.4) were discussed in detail in [10, 11, 13]. The 
results of the calculations in the framework of the IVM the leading-order are presented in 
Appendix A. 

Since both <J>n (a, q 2 = ml) « <&m (a, 0) and 4>2;n(&) are normalized to unity, we im- 
mediately conclude that 

Ia = >/2fl (4.9) 

Indeed, straightforward calculation shows that the relation (4.9) is satisfied with accuracy 
of about ~ 5%, which is mainly related to the precision of numerical evaluation. 

In Figure 10 we present the ratio of the left and right hand sides of Eqs. (4.5) - (4.8) 
in order to demonstrate how accurate is PCAC at different values of a. As we can see, 
PCAC is valid for all functions. The ratio 4>2\tv/^\\ near the endpoints a — > 0, 1, vanishes 
because both functions (/>2 ;7r and $n are suppressed at the endpoints, and as a result we 
have large numerical errors in those regions. For the ratio gs/ipA-^ both DAs have very 
similar shapes. However, the position does not exactly match: numerical evaluation yields 
d^ 4 . v ~ 0.786, ctg 3 ~ 0.791. For all the other functions the PCAC relation is exact. 

5 Conclusion 

Partial conservation of axial current is a fundamental hypothesis, which involves the non- 
trivial dynamics of spontaneous breaking of chiral symmetry. This hypothesis has been 
tested in various process at low energy [3] , and in the process of diffractive neutrino inter- 
action at high energies [24, 25, 29]. We perform a direct test of the PCAC relation between 
the derivative of the axial current and the pion field, by comparing the distribution ampli- 
tudes, calculated within the light-cone dipole approach. Because PCAC essentially involves 
soft interactions, a perturbative QCD technique, usually used for calculation of the distri- 
bution amplitudes, cannot be employed. We derived several relations between the axial 
current and pion distribution amplitudes, and performed calculations in the framework of 
the instanton vacuum model. Our results confirmed that the based on PCAC relations are 
correct with high accuracy. 
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Figure 10. Left: Numerical verification of the PCAC relations within the IVM. According to 
PCAC, all the ratios should be equal to one. Right: comparison of g% and "04 ;7r - Because gs and 
■04;7r have zeros at the points a 93 ~ 0.5 ±0.291 and a^, 4 . T « 0.5 ±0.286, the ratio plotted by dashed 
line in the left figure is significantly different from one near these points. 



A Pion distribution amplitudes 

The distribution amplitudes of the pion (4.2)-(4.4) and their evaluation in the framework 
of the IVM was performed and discussed in detail in [10, 11, 13, 15-17]. For the sake 
of completeness we present here the leading-order expressions for the pion distribution 
amplitudes. 
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